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LIST  OF  SYMBOLS  i 


a constant  in  Tillotson’s  equation  of  state 
A constant  in  Tillotson's  equation  of  state 
b constant  in  Tillotson's  equation  of  state 
B constant  in  Tillotson's  equation  of  state 
c sound  speed 

c^  constant  in  sound  speed  equation 
Cj^  constant  in  sound  speed  equation 


d°  penetration  associated  with  t° 

f f 

dj  penetration  associated  with  tj 

d^  penetration  associated  with  t^ 

G modulus  of  rigidity 

I specific  internal  energy 

Ij^  constant  in  yield  strength  equation 

constant  in  Tillotson's  equation  of  state 

Ig  constant  in  Tillotson's  equation  of  state 

constant  in  Tillotson's  equation  of  state 


decrease  in  rod  or  jet  length 

steady-state  ratio  of  rod  loss  to  penetration  in  the  plate 

steady-state  ratio  of  rod  loss  to  penetration  in  the 
semi-infinite  material 


p pressure 

p^  pressure,  condensed  state 

Pg  pressure,  expanded  hot  gases 

P depth  of  penetration  or  target  thickness 

r radial  coordinate 


p 


t time  for  rarefaction  wave  to  travel  the  radius  of  the  rod 
0 

tj^  time  to  penetrate  12.7  mm  (^s  in.)  of  plate  material 

t^  time  to  penetrate  25.4  mm  (1  in.)  of  plate  material 

tj  time  at  the  beginning  of  Vj 

t°  tim-r  at  the  beginning  of  V2 

f 

tj  time  at  the  end  of 

impact  velocity 
Vj  jet  velocity 

Vp  . penetration  velocity 
Vj^  steady-state  velocity  in  the  plate 

V2  steady-state  velocity  in  the  semi-infinite  material 

V yield  strength  in  shear 

constant  in  yield  strength  equation 

Yj  constant  in  yield  strength  equation 

Y2  constant  in  yield  strength  equation 

z longitudinal  coordinate 

a constant  in  Tillotson's  equation  of  state 

5 constant  in  Tillotson's  equation  of  state 

Ar  radial  grid  increment 

Az  longitudinal  grid  increment 

V n - 1 

V consta.it  in  the  sound  speed  equation 

p mass  density 

pj  mass  density  of  the  jet 

p ambient  mass  density 

p^  mass  density  of  the  target 
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I . INTRODUCTION 


This  report  is  concerned  with  the  behavior  of  a shaped-charge  jet 
when  it  penetrates  a target  possessing  a material  discontinuity,  for 
example,  a steel  cover-plate  on  top  of  a thick  aluminum  plate.  This 
study  is  limited  to  targets  whose  cover-plate  material  is  more  dense 
than  the  second  material.  A material  discontinuity  occurs  at  the 
interface  between  the  two  materials,  and  we  shall  refer  to  any 
interaction  of  the  jet  with  this  interface  simply  as  rod-target 
interaction.  A computer-oriented  mathematical  model  was  used  in  this 
study  of  rod-target  interaction. 

There  are  several  two-dimensional  models  (codes)  that  were  developed 
to  handle  penetration  problems  and  are  operational  at  the  Ballistic 
Research  Laboratory.  These  codes  are  OIL,^  TOIL,^  DORF,^  and  HELP,**»^ 
which  are  Eulerian  formulated,  and  HEMP,^  which  is  Lagrangian  formulated. 
The  large  deformation  of  material  occurring  during  a shaped-charge 
penetration  influenced  the  selection  of  one  of  the  Eulerian  codes  for 
this  study,  and  of  these,  the  HELP  code  was  selected  because  of  its 
capability  of  denoting  the  material  interfaces  and  free  surfaces. 

The  objective  of  this  computer  study  is  to  investigate  the 
parameters  affecting  rod-target  interaction  in  order  to  provide  data  on 
penetration  and  rod  loss  that  could  be  used  by  armor  experimentalists 
and  designers.  The  results  of  the  computer  study,  which  are  reported 
here,  include  only  that  part  of  the  voluminous  data  generated  by  HELP 
that  pertains  to  penetration  depth  and  loss  of  rod  material. 


/.  W.  E.  JohnAon,  ”01 L,  A ContimouA  T(M~V-bnznAiomZ  EuleAMin 

Hyd/iodynmic.  Code,”  Gul^  GzncAal  Atomic,  GAMP-5580  Reused,  JamoA.y 
1965. 
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W.  E.  JohnAon,  "TOIL,  A Tux)~Mate/UaZ  VeAAion  the  OIL  Code,”  GuZ{ 
GeneAot  Atomic,  6AMP-S073  Addendum,  November.  1967. 

WaJtla.ee  E.  JohnAon,  "Code  ConAeZxvUon  Study,"  Ain.  FoA.ee  WeaponA 
LaboAjotony  Technical  RepoAt  Uo.  70-144,  Apiii  1971. 

L.  J.  Hageman  and  J.  M.  WatAh,  "HELP,  A Multi-MateAial  Eule/Uan 
PAogAom  (JoA  CompA.eAAible  Fluid  and  EloAtic-PlaAtic  FIoica  in  Two 
Space  VimenAionA  and  Time,  l/olume  I,"  BalliAtic  ReAeoAch 
LabonatoAieA  ContAact  RepoAt  No.  39,  May  1971.  (AD  #726459) 


L.  J.  Hageman  and  J.  M,  WalAh,  "HELP,  A Multi- MatcAial  EuleAian 
PAogAom  iJoA  CompACAAible  Fluid  and  EloAtic-PloAtic  FIowa  in  Two 
Space  VimenAioM  and  Time,  Volume  II:  FORTRAN  LiAting  o^  HELP," 
BalliAtic  ReAea/ich  LaboAotoAieA  ContAact  RepoAt  No.  39,  May  1971. 

"CalculationA  OjJ  EloAtic-PloAtic  Flow,"  UniveAAity 
0(J  Cali^oAnia,  Lamence  Radiation  LaboAotoAiy,  UCRL-7322,  ApAil  1963. 
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II.  COMPUTER  CODE 
A.  Description  of  the  HELP  Code 

HELP  Can  acronym  for  Hydrodynamic  Elastic  Plastic)  is  a FORTRAN 
code  developed  by  Systems,  Science  and  Software,  Inc.  It  evolved  from 
three  major  hydrodynamic  codes  previously  developed  over  a 20-year 
period.  In  chronological  order  of  development,  these  codes  are  PIC,^ 

OIL,  and  RPM.® 

This  code  is  Eulerian  formulated,  that  is,  a numerical  mesh  or  grid 
is  fixed  in  space  and  material  is  permitted  to  flow  or  move  through  the 
grid.  Material  motion  is  governed  by  the  numerical  approximations  to 
the  equations  of  conservation  of  mass,  momentum,  and  energy.  These 
equations  along  with  an  equation  of  state  and  a stress-strain  relation- 
ship provide  sufficient  equations  for  their  numerical  solution. 

HELP  is  time  dependent  in  two  space  dimensions  with  an  option  for 
either  plane  (x,y)  coordinates  or  cylindrical  (r,z)  coordinates.  When 
the  cylindrical  coordinate  option  is  used,  the  problems  to  be  solved 
necessarily  involve  axisymmetric  configurations  such  as  spheres,  cones, 
infinite  plates,  and  right  circular  cylinders  which  impact  normally 
upon  axisy?nmetric  targets.  The  time  increment,  which  is  determined  by 
the  code,  is  based  on  the  Courant  stability  condition  and  the  maximum 
sound  speed  plus  the  particle  velocity  in  the  grid. 

The  number  of  materials,  which  can  be  treated  in  one  problem,  is 
limited  only  by  the  size  of  the  memory  of  the  computer  being  used.  The 
BRL  versior  of  HELP  can  handle  three  materials  in  addition  to  a vacuum. 
The  different  materials  are  separated  from  each  other  by  lines  joining 
massless  tracer  particles.  These  particles  move  across  the  grid  and 
locate  the  interfaces  and  free  surfaces  of  each  material.  Consequently, 
a Lagrangian  effect  is  introduced  into  the  treatment  of  moving  surfaces. 

The  HELP  code,  using  the  von  Mises  yield  condition,  treats 
materials  as  being  elastic-plastic. 

The  input  for  HELP  involves  the  specification  of  initial  density, 
specific  internal  energy,  and  velocity  in  each  cell  containing  any 
material.  Furthermore,  locations  of  tracer  particles  along  free  surfaces 
and  interfaces  are  specified. 


7.  Anthony  A.  Amdtn,  "The  PoAticZe^tn^Cett  Method  ^OA  the  VyncuiUcs  o^ 
CompAeiitbZe  Fhud6,"  Loa  Alcmo6  Sctenti^tc  LaboAotoAy  oi  the 
UntoeAJiity  oi  CaJLiloAtuA,  LA-3466,  PebAuaAy  1966, 

t,  J.  K,  Vtenei,  M,  U.  EvanA,  L,  J.  Hageman,  W,  Z,  Johnson,  and  J,  M. 
(ilalsh,  "An  EuZeAZan  Method  {oA  CaZeuZating  StAength  Dependent 
Ve^oAmatZon,  PoAt  One,  PoAt  Two,  PoAt  ThAee,  and  Addendum,"  GuZi 
GeneAoZ  Atomic,  GAMD-S497,  Fd^AuoAy  1968. 
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B.  Computer  Output 


The  output  of  the  BRL  version  of  HELP  lists,  at  specified  time 
intervals,  the  total  internal  energy,  kinetic  energy,  mass,  momentum 
components,  and  plastic  work  done  for  each  material.  The  total  energy, 
calculated  for  each  computational  cycle,  must  agree,  within  specified 
limits,  with  the  total  theoretical  energy,  otherwise,  the  computational 
process  is  brought  to  a halt. 

The  coordinates  of  the  tracer  particles  are  listed  as  well  as,  for 
each  cell  containing  any  material,  the  velocity  components,  hydrostatic 
pressure,  specific  internal  energy,  density,  mass,  and  stress  components. 


Plotting  programs  have  been  developed  at  BRL  to  display  the 
voluminous  data  generated  during  a typical  computer  run.  Examples  of 
the  types  of  plots  which  can  be  generated  are  shown  in  Figures  1,  2,  and 
3.  These  plots  were  generated  from  data  computed  in  this  study; 
however,  this  type  of  information  will  not  be  reported  here  except  for 
a brief  description  of  these  figures. 

Figure  1 shows  the  pressure  field  at  a particular  time  after  the 
impact  of  a right  circular  rod  on  a plate  backed  by  a semi-infinite 
material;  Figure  2,  the  velocity  field;  and  Figure  3,  the  stress  field. 
In  each  of  these  figures  the  tracer  particles  were  used  to  outline  the 
interfaces  and  free  surfaces  of  the  different  materials.  Since  the 
problem  is  axisymmetric,  only  half  of  the  various  fields  are  presented. 

C.  Equation  of  State 

Since  HELP  has  an  equation-of-state  subroutine,  one  may  incorporate 
any  equation  of  state  into  the  program  that  is  of  the  form  where  the 
pressure  is  a function  of  mass  density  and  specific  internal  energy. 

The  BRL  version  of  HELP  contains  an  equation  of  state  for  an  ideal  gas 
and  also  Tillotson's  equation  of  state. ^ 

Tillotson's  equation  of  state  takes  two  forms: 

1.  For  the  condensed  states,  that  is,  when  p > p , or  for 
any  cold  states,  that  is,  when  I < I^, 


p = Pc  = 


a ♦ 


♦ 1 


Ip  Ap  ♦ Bp‘ 


2.  For  expanded  hot  gases,  that  is,  when  P and  I > I^, 


p = p = alp  + 


bip 


1 


'o'’ 


+ Ape 


-B(l/n  - 1) 


-a(l/n  - D' 


9. 


J.  H.  TlUotion,  "MeXatUc  Equations  State,  ion  HypeAveZocAjty 
Impact,"  Guti  GeneAol  Atomic,  GA-3216,  July  1962. 
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where 

p = pressure 

I « specific  internal  energy 
p s mass  density 
n * p/pQ 

u » n - 1 

and  a,  b.  A,  B,  a,  B,  p^,  I^,  I^,  and  are  constants  for  a particular 

material.  The  values  of  these  constants  for  materials  which  were  used 
in  this  report  are  listed  in  Table  I. 

In  addition  to  Tillotson's  equations,  HELP  provides  a transition 
equation  between  the  condensed  and  expanded  states  where  p < p and 

T ^ T ^ T I . O 


Cl  - I3)  Pe  ^ (V  - I)  P^ 


Pressure  versus  particle-velocity  Hugoniots  and  shock  velocity 
versus  particle  velocity  graphs,  based  on  Tillotson's  equations,  are 
shown  in  Figures  4 and  5,  respectively. 

D.  Strength 

Several  more  material  constants  are  required  as  inputs  for  HELP. 
The  modulus  of  rigidity,  G,  relates  the  deviatoric  stress  to  the 
deviatoric  strain  by  a linear  elastic  relation.  A variable  yield 
strength  in  shear,  defined  by 

y . (Y„  * YjU  . YjW^)  Cl  - i/y, 

requires  four  additional  constants.  The  values  of  these  constants,  for 
the  materials  used  in  this  report,  are  given  in  Table 


III.  ROD-TARGET  PARAMETERS 

The  configuration  of  the  jet  and  target  for  the  computer  study  of 
rod-target  interaction  consisted  of  a copper  rod,  or  jet,  impacting 
normally  upon  a steel  plate  which,  in  turn,  was  backed-up  by  a semi- 
infini  e material.  The  parameters  for  this  study  are  listed  in  Table 
III. 

The  length  of  the  rod  was  chosen  to  be  long  enough  to  prevent 
rod-end  effects,  such  as  shock  reflections  or  rarefactions,  from 
affecting  the  rod-target  interaction.  For  simplicity,  the  rod,  prior  to 
impact,  was  given  a uniform  axial  velocity;  actually,  a velocity  gradient 
exists  in  a shaped  charge  jet. 

10.  William  Giltich  and  GeoAge  E.  HauveA,  JnleA.-o^(ice.  coAAMpondence, 
FcbfuuULtf  1970, 
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It  was  not  intended  that  all  combinations  of  the  parameters  listed 
in  Table  III  would  be  run  on  the  computer;  rather,  the  selection  of 
combinations  of  parameters  would  be  guided  by  the  results  of  sets  of 
computer  runs.  The  runs  were  made  at  the  Kirtland  Air  Force  Base  on 
the  CDC  6600.  The  matrix  of  parameters  for  these  runs  is  listed  in 
Table  IV  in  the  order  in  which  they  were  run. 


IV.  THEORY  OF  PENETRATION 

Before  presenting  the  numerical  results  obtained  from  the  computer 
runs,  we  shall  present  some  preliminary  calculations  of  the  penetration 
process  for  use  as  a comparison  with  the  numerical  results. 

When  a copper  rod  impacts  upon  a steel  plate,  shock  waves 
immediately  propagate  from  the  interface  into  the  steel  plate  and  the 
copper  rod.  For  a short  time  after  impact,  the  high  pressure  in  the 
region  along  and  near  the  axis  of  symmetry  of  the  rod-plate  configuration 
and  between  the  two  diverging  shocks  remains  constant.  Eventually, 
rarefactions  from  the  front  surface  of  the  plate  and  the  sides  of  the 
rod  will  relieve  this  pressure.  Estimates  of  the  initial  impact 
pressure  and  particle  velocity  can  be  obtained,  by  the  method  of  reflected 
Hugoniots,^^  from  Figure  4.  Thereafter,  Figure  5 can  be  used  to 
determine  the  shock  velocities.  A summary  of  these  values  is  presented 
in  Table  V. 

The  sound  speed,  c (km/s),  in  the  high  pressure  region  is  determined 
from  the  relation^ ^ 

c = Cjj  + Cjp^ 

where  the  constants  c^  (km/s),  c^,  and  v are  3,20,  3.405,  and  0.452, 

respectively,  for  copper  and  3.50,  3.259,  and  0.439,  respectively,  for 
iron.  The  pressure,  p,  is  in  units  of  megabars. 

An  estimate  of  the  time,  t^,  required  for  the  lateral  rarefactions, 

coming  from  the  edge  of  the  contact  surface,  to  reach  the  axis  of 
symmetry  of  the  rod-plate  configuration  is  made  by  dividing  the  radius 
of  the  rod  by  the  sound  speed.  Table  VI  sumarizes  this  information. 

Since  the  sound  speed  is  greater  in  iron  than  copper,  t was  calculated 
from  the  sound  speed  of  iron.  ° 

In  the  computer  runs  that  were  made,  the  127-ram  (5-in.)  length  of 
rod  and  the  12.7-  mm  (*s-in.)  and  25,4-mm  (1-in.)  plate  thicknesses  were 
large  enough  to  prevent  any  disturbances,  from  the  back  of  the  plate  or 

11.  R.  CouAant  and  K.  0.  Tfu.tdAA.chi>,  Sup^otUc  Tim  and  Shock  Uavti, 
Volumt  I,  New  VoAk,  JnttAAcitnet  PubCcihtu,  Inc.,  1948. 

IS.  ChtittA  GA06ch,  "Sound  VcZocaXaca  Along  the  Hugoniot  ioK  Hint 
Melals  cU  High  Pacmuac,"  Sttvtni  iMtUatt  0(J  Ttchnology, 

UnpubtLshtd  note,  JanuoAy  196S. 


the  rear  of  the  rod,  from  reducing  the  duration  of  constant  pressure 
listed  in  Table  VI. 

Following  the  constant  pressure  phase  of  the  penetration  process, 
a steady-state  phase  begins.  Here,  the  rate  of  penetration  (the 
velocity  of  the  rod-plate  interface  at  the  axis  of  symmetry)  remains 
constant.  This  phase  continues  until  disturbances  from  the  interface 
of  the  plate  and  semi- infinite  material  interact  with  the  main  shock 
in  the  target. 

A transient  phase  follows  during  which  the  rod- target  interaction 
occurs.  This  phase  is  followed  by  another  steady-state  phase  in  which 
the  penetration  rate  has  adjusted  itself  to  the  semi-infinite  material. 

Estimates  of  the  penetration  velocity  and  the  decrease  in  rod 
length  during  the  steady-state  phases  can  be  made  from  the  following 
relations  based  on  Bernoulli's  equation  of  incompressible  flow:^^ 


1 ♦ (P^/Pj) 


2.  P - L(p 


where 


Vp  = penetration  velocity 
V,  * jet  velocity  at  impact 

p^  s density  of  target 
Pj  = density  of  jet 

P = depth  of  penetration  or  target  thickness 
L - decrease  in  jet  length 

Penetration  velocities  were  calculated  from  Equation  1 for  three 
jet  velocities  where  the  jet  material  was  copper  and  the  target  material 
was  varied.  Assuming  that  the  penetration  velocity  remains  constant  as 
the  jet  penetrates  a steel  target,  the  times,  tj^  and  tj,  required  to 

penetrate  12.7  mm  (*s  in.)  and  25.4  mm  (1  in.)  of  steel,  respectively, 
were  also  calculated.  These  data  are  listed  in  Table  VII. 

Figures  6,  7,  and  8 show  the  displacement  along  the  axis  of 
symmetry,  as  a function  of  time,  of  the  jet-target  interface  at  impact 
velocities  of  3,  5,  and  7 km/s,  respectively.  Each  target  was  a 
25.4-mm  (1-in.)  plate  of  steel  backed  with  various  materials.  These 
graphs  were  generated  from  the  information  in  Table  VII  under  the 
assumption  that  the  rod  was  long  enough  to  prevent  end  effects.  Similar 
graphs  could  be  plotted  for  a 12.7-rora  (4-in.)  plate. 
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V.  NUMERICAL  RESULTS 

Although  a vast  c,  intity  of  data  is  generated  with  each  computer 
run  of  the  HELP  code,  only  data  on  the  motion  of  the  free  surface  and 
interfaces  are  reported  here,. 

For  every  cycle  of  'omputation,  the  positions  of  three  tracer 
particle.",  were  recorded.  These  tracer  particles  were  located  along  the 
axis  of  symmetry  at  the  oa.se  of  the  rod,  the  rod-plate  interface,  and 
the  'Interf-ice  of  the  pif.:.e  and  semi-infinite  material.  Figure  9,  which 
depicts  the  grid  and  pr,»;(e<.J  ile-tar;je-t  configuration  used  in  these 
computer  runs,  shows  t’,ie  ;.->cations  of  these  tracer  particles. 

»iith  reference  to  tie  grid  i '^’gure  9,  the  r-increment  (the  length 
of  each  cell  in  the  r-diroctio.i) , .ir,  is  equal  to  0.0079375  mm  for  the 
first  sixteen  colls  or  c ^u.  ns,  ’ hei  lafter,  the  r-increment  was  increased 
geometrically,  with  a ratio  j.i,0.roxln.,\tely  1.1,  until  a total  of  forty 
columns  wore  dimensioned  in  the  r-direction.  The  developer  of  HELP 
recommends  that  the  geometric  ratio  does  not  exceed  1.5,  The  r-increment 
w?.s  i.'iv  e.?sed  in  the  r-,_>ection  to  piovide  space  for  the  shock  waves, 
'.'"seating  t' the  iuptu. ..  region,  zo  travel  and  be  analyzed  before 
reaching  'ne  riglit  boundary  of  the  t‘f»**’itih  '.Dl.i.-n,  This  boundary  is 
tran.Mnis  ;ive,  that  is,  *Ti.  js  ma^‘  flow  xhr.ugh  this  boundary  and  out  of 
t.’.K-.  gri.,. 

T.ie  g" -a  consist-,  of  ells  or  rows  in  the  z-direction  with  the 
/-increment  of  thy  tc>.  8S  lu  » being  equal  to  0,018143  mm;  thereafter, 
the  ,:‘  i;icr('m&iit,  Az,  was  inc- cased  geometrically  over  14  rows  iu  the 
negative  z-direction  with  s lario  of  approximately  1,1376,  This  increase 
in  the  z-increment  provided  sp&cv;  cor  the  long  rod.  The  bottom  row 
had  a Az  of  0.018143  mm;  however,  this  dimension  is  not  important  since 
there  was  no  material  motion  in  this  row.  The  top  and  the  bottom  of 
the  grid  were  transm.'isive  diso. 

The  aspect  ratio,  Lz/hr,  was  approximately  2,286  in  the  region 
where  the  r and  z increments  were  constant.  We  would  have  preferred  that 
the  aspect  ratio  did  not  exceed  1..5,  as  recommended  by  the  developer  of 
HELP;  however,  the  long,  thin  rod  and  the  thickness  of  tb**  target 
r-ohioited  the  use  of  a smallrr  aspect  ratio.  The  radius  of  the  rod 
w.-"*.  ’’epresented  by  only  fou-'  cells  which  results  in  poor  resolution  of 
S'v  ek  and  rarefaction  fronts.  If  a larger  number  of  columns  were  used 
fer  the  rod's  radius,  the  Ar  would  have  to  be  smaller  and  the  aspect 
ratio  would  be  increased. 

With  the  constant  Az  that  was  used,  the  thickness  of  the  25.4-iilm 
(l-in.)  plate  was  represented  by  14  rows  of  cells  and  by  7 rows  of  cells 
for  the  12.7-mm  (^s-in.)  plate.  Since  a shock  front  is  usually  smeared 
over  4 or  ' c.-’lls,  we  should  not  expect  good  resolution  of  the  shock  front 
in  th"  z-direction.  The  number  of  rows  representing  the  plate  could  have 
He*  increased  by  increasing  the  number  of  rows  in  the  grid;  however,  this 
would  increase  the  running  time  and  computer  memory  requirements  for  the 
the  problem.  With  the  grid  that  was  used,  the  running  time  for  each 
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case  for  approximately  three  hours  with  the  use  of  81,000  words  of 
memory,  which  was  near  the  capacity  of  the  CDC  6600,  Consequently,  with 
the  grid  that  was  used,  sharp  shock  detail  was  not  expected;  rather,  a 
general  picture  of  the  penetration  phenomena  was  expected. 

The  displacement,  as  a function  of  time,  of  the  ti'acer  particles, 
shown  in  Figure  9,  were  graphed  in  Figures  10-21  for  the  cases  that  were 
run  on  the  computer.  This  information  is  not  shown  for  Case  10  as  there 
was  difficulty  in  reading  the  data  from  the  magnetic  tape  output  of  the 
computer.  These  graphs  represent  the  prime  output  of  the  runs  that  were 
made.  The  following  sections  provide  a discussion  of  these  graphs  and 
further  information  obtained  from  them. 


VI.  DISCUSSION 

For  all  the  cases  that  were  run  on  the  computer,  the  base  of 
the  rod  moved  at  a constant  velocity  over  the  duration  of  the  run  as 
evidenced  by  the  constant  slope  of  the  trajectory  of  the  base  tracer 
particle  in  Figures  10-21.  This  indicates  that  either  none  of  the 
disturbances,  initiated  at  the  front  end  of  the  rod,  reached  the  base  of 
the  rod  during  the  time  period  considered  or  the  code  could  not  detect 
these  disturbances. 

The  interface  of  the  plate  and  the  semi-infinite  material  begins  to 
move  before  the  rod-plate  interface  reaches  it  since  a shock  wave 
precedes  the  moving  interface.  Note  that  the  thickness  of  the  steel 
plate  along  the  axis  of  symmetry  decreases  with  time,  but  never  becomes 
zero,  that  is,  the  plate  is  not  perforated.  This  derives  from  the 
fact  that,  with  the  HELP  code,  a tracer  particle,  initially  positioned 
on  the  axis  of  symmetry,  cannot  move  off  the  axis,  but  only  along 
the  axis.  Furthermore,  these  tracer  particles  cannot  pass  one  another. 
Thus,  some  plate  material  will  always  separate  the  rod  material  from 
the  semi-infinite  material.  In  actuality,  the  plate  would  be  perforated. 

The  displacement-time  data  from  the  computer  output  was  interpolated 
to  obtain  tj^  and  t^,  the  time  expended  while  the  rod-plate  interface 

tracer  particle  moved  12.7  mm  (b  in.)  and  25.4  mm  (1  in.),  respectively 
(depending  on  the  thickness  of  the  plate).  See  Table  IX  for  a tabulation 
of  these  times.  For  a given  velocity,  these  times  increase  as  the  density 
of  the  semi-infinite  material  is  increased,  A comparison  of  the  data  for 
an  all  steel  target  (where  the  semi-infinite  material  is  steel)  with  the 
data  in  Table  VII  indicates  that  t,  and  t,,  obtained  from  Bernoulli's 

incompressible  equation,  are  less  than  that  obtained  from  the 
compressible  material  code,  HELP. 

If  we  measure  penetration  from  the  initial  position  of  the  free 
surface  of  the  steel  plate,  the  displacement  of  the  rod-plate  tracer 
particle  gives  us  penetration.  Penetration  trajectories  are  shown  in 
Figures  22  and  23  for  plate  thicknesses  of  12.7  mm  (3s  in.)  and  25.4  mm 
(1  in.),  respectively.  One  may  wish  to  compare  Figure  23  with  Figures  6, 
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7,  and  8 to  see  the  difference  between  the  compressible  and  incompressible 
mathematical  models  of  penetration. 

Examining  Figure  23,  we  see  that  the  penetration  process  starts 
with  a fluctuation,  the  duration  of  which  appears  to  decrease  as  the 
impact  velocity  is  increased.  The  duration  of  this  period  is  denoted 

by  t°.  A steady-state  period,  characterized  by  a constant  slope 

(velocity),  of  the  penetration  trajectories,  follows  and  lasts  until 

f 

a time  denoted  by  tj^. 

f 

Following  tj,  the  penetration  trajectories  pass  through  a 
transition  period  ending  at  time  t^,  which  marks  the  beginning  of 
another  steady-state  period.  Rod-target  interaction  occurs  between 
tj  and  t°;  this  will  be  discussed  later. 

To  gain  further  information  on  the  penetration  process,  let  us 
examine  the  penetration  velocity  histories  of  the  cases  which  were  run 
on  the  computer.  The  penetration  velocity  - time  curves,  shown  in 
Figures  24  and  25,  were  obtained  by  fitting  cubic  spline  functions^**  to 
the  data  in  Figures  22  and  23  and  taking  the  first  derivatives  of  the 
functions.  Using  the  data  in  Figures  22  and  23,  we  have  also  plotted 
penetration  velocity  - penetration  curves  as  shown  in  Figures  26  and 
27.  These  velocity  curves  show  that  the  penetration  process  starts  out 
with  a fluctuation  in  velocity,  as  pointed  out  previously. 

The  steady-state  velocity,  Vj,  between  t°  and  tj,  was  determined 

by  a least-squares  fit  to  the  data  that  was  close  to  being  a straight 
line.  This  velocity  is  listed  in  Table  IX.  The  steady-state  velocity, 

V2,  which  occurs  after  t^,  was  similarly  determined  and  is  also  listed 

in  Table  IX.  The  values  for  and  ^ tend  to  average  about  6%  less 

than  the  values  listed  for  V in  Table  VII. 

P 

The  values  of  the  times,  t°,  t°,  and  t^,  which  are  tabulated  in 
Table  IX,  were  determined  by  estimating  when  Vj  and  V2  began  and  when 

Vj  ended.  If  we  compare  t^  in  Table  VI  with  t°  in  Table  IX,  we  find 

that  it  takes  an  order  of  magnitude  more  time  to  reach  a steady-state 
velocity  than  for  a rarefaction  wave  to  reach  the  axis  of  symmetry  of 
the  rod-target  configuration. 


14,  PaJbmA  R.  SMegeZ,  "AppxoxAjnatz  by  Cubic  Sptinc  uoijth  Rcip&ct  to 
Eactidcan  Uohm/*  BaZtUtic  ReieoAch  LuboKcvtonJizA  Rcpofit  No.  1465, 
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Figures  22  and  23  were  used  to  determine  d 
penetration  distances  corresponding  to  times,  t 


These  distances  are  listed 


1*  ^ 
in  Table 


and  d!,,  the 


1» 

IX. 


and  t 


2‘ 


respectively. 

The  rod  loss,  L,  was  determined  by  subtracting  the  distance  between 
the  base  tracer  particle  and  the  rod-plate  interface  tracer  particle 
(Figures  10-22)  from  the  initial  length  of  the  rod,  127  mm  (5  in.).  This 
information  is  presented  as  rod  loss  - penetration  curves  in  Figures 
28-31.  By  a least-squares  fit  of  the  data  in  these  figures  that  were 
close  to  a straight  line,  we  determined  the  rod  loss  to  penetration 

ratios  for  the  plate  material  and  the  semi-infinite  material,  and 


respectively.  These  ratios,  which  are  listed  in  Table  IX,  are 


greater  than  the  values  of  ^ in  Table  VIII  for  incompressible  flow. 


Between  the  penetration  distances  of  d^  and  d°,  the  slope  of  the 

curves  in  Figures  28-31  is  changing.  This  is  the  transition  region 
where  the  penetration  process  is  affected  by  the  unlike  material 

£ 

interface  in  the  target.  Between  dj  and  a penetration  distance  equal 

to  the  original  thickness  of  the  plate,  the  gradient  of  the  rod  loss, 
with  respect  to  penetration,  starts  to  decrease.  Guided  by  Equation  2, 
we  could  say  that  the  plate  material  was  acting  as  if  it  were  less  dense 
than  normal.  After  a penetration  distance  greater  than  the  original 

plate  thickness,  and  less  than  d^,  is  reached,  the  rod-loss  gradient, 

with  respect  to  penetration,  is  greater  than  that  finally  obtained  in 
the  steady-state  region  in  the  semi-infinite  material.  Again  according 
to  Equation  2,  we  could  say  that  the  semi-infinite  material  was  acting 
as  if  it  were  more  dense  than  normal. 


From  Figures  10-22,  we  observe  that,  not  only  is  the  thickness  of 
the  steel  plate  (the  distance  between  the  rod-plate  tracer  particle  and 
the  tracer  particle  for  the  interface  of  the  plate  and  the  semi-infinite 
material)  decreasing  along  the  axis  of  symmetry,  but,  more  important, 
some  of  the  p’ate  material  has  been  pushed  into  the  region  initially 
occupied  by  the  semi-infinite  material.  For  example,  Figure  10  indicates 
that  approximately  4.06  mm  (0.16  in.)  of  plate  material  has  moved  into 
this  region  when  the  rod  has  penetrated  a distance  of  25.4  mm  (1  in.), 
the  original  thickness  of  the  plate.  This  is  one  of  the  reasons  that  the 
semi-infinite  material  is  acting  as  if  it  were  more  dense  than  normal. 


A method  that  we  could  use  to  gauge  the  relative  merits  of  various 
materials  in  hindering  penetration  is  by  using  Equation  2 to  determine 
an  apparent  target  density.  The  density  of  the  rod  would  be  the  initial 

copper  density  and  the  ratio,  would  be  determined  from  the  numerical 
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results.  However,  since  the  penetration  velocity  histories  have  already 
been  determined,  Equation  1 will  be  used  for  this  purpose.  This 
analysis  is  shown  in  Figures  32-35,  where  the  apparent  density  has  been 
normalized  with  respect  to  the  initial  density  of  steel  for  the  first 
25.4  mm  (1  in.)  (the  initial  thickness  of  the  plate)  of  penetration; 
thereafter,  the  target  density,  obtained  from  Equation  1,  has  been 
normalized  with  respect  to  the  initial  density  of  the  semi-infinite 
material.  Of  course,  this  normalization  is  undefined  for  a vacuum. 

Examining  Figure  33,  we  again  see  fluctuations  in  the  curves  during 
the  initial  penetration  process.  Next  comes  the  steady-state  region 
where,  as  we  previously  pointed  out,  the  penetration  velocity  obtained 
from  the  numerical  results  was  less  than  that  determined  from  the 
incompressible  flow  theory.  This  difference  in  velocity  accounts  for 
the  normalized  density  being  greater  than  unity  as  one  would  expect. 
Nevertheless,  some  information  on  the  effect  of  the  material  of  the 
semi-infinite  medium  on  rod  loss  can  be  obtained  from  this  analysis 
when  we  look  ' transition  region  between  the  two  steady-state 

regions. 

Just  t ;-f;,re  25.4  mm  (1  in.)  of  penetration,  we  see  that  the 
semi-infinite  material  with  the  lower  initial  density  (aluminiun  compared 
to  titanium)  provides  a lower  apparent  density  for  the  steel  plate. 

Thus,  the  less  dense  semi-infinite  medium  causes  less  rod  loss  in  this 
region  than  a more  dense  medium.  However,  just  after  25.4  mm  (1  in.)  of 
penetration,  the  semi-infinite  material  with  the  lower  initial  density 
(aluminum  compared  to  titanium)  provides  a higher  apparent  density  for 
the  semi-infinite  medium.  Thus,  the  less  dense  semi-infinite  medium 
causes  more  rod  loss  in  this  region  compared  to  a more  dense  medium. 

An  apparent  density  curve  for  air  as  a semi-infinite  material  was 
not  plotted  since  its  low  density,  compared  to  the  density  of  the  metals 
that  were  used  as  semi-infinite  material,  resulted  in  apparent  densities 
that  were  three  orders  of  magnitude  greater  than  shown  in  Figures  32-35. 

A comparison  of  the  steel-aluminum  curves  in  Figures  32-34  does 
not  indicate  any  clear  effect  of  impact  velocity  on  tht  apparent  density 
of  the  semi-infinite  material. 


VII.  SUMMARY 

The  results  from  a set  of  computer  runs,  designed  to  investigate 
the  interaction  of  a copper  rod,  simulating  a shaped-charge  jet,  with 
the  interface  of  the  target  materials,  were  presented.  The  results  were 
mainly  concerned  with  the  penetration  process  as  viewed  along  the  axis 
of  symmetry  of  the  rod-target  configuration. 

A comparison  of  the  numerical  (unsteady,  compressible)  results  with 
the  results  from  Bernoulli’s  incompressible  equations  shows  that  the 
latter  results  may  be  satisfactory  for  predicting  rod  loss,  penetration 
velocity,  etc.  for  many  purposes;  however,  for  an  experimentalist  or 
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armor  designer  to  take  full  advantage  of  phenomena  which  can  hinder 
P'^netration,  compressibility  should  be  taken  into  account. 

The  numerical  results  show  that  the  penetration  process  starts  with 
a fluctuation  in  the  penetration  velocity,  followed  by  a steady-state 
period.  Naturally,  this  steady-state  period  will  occur  only  if  the 
target  plate  is  thick  enough.  The  results  fr  m the  runs  of  the  12.7-mm 
C^-in.)  plate  indicate  that  this  thickness  was  Just  about  on  the 
borderline  for  not  getting  a steady-state  condition.  Following  the 
steady-state  period,  a transition  period  occurs  as  the  penetration 
velocity  adjusts  to  the  final  steady-state  velocity. 

This  transition  region  can  be  divided  into  two  zones;  one  before  a 
penetration  distance,  equal  to  the  initial  thickness  of  the  plate,  is 
reached  and  one  after  this  depth  is  reached. 

In  the  first  zone,  the  apparent  density  of  the  steel  is  lower  than 
its  normal  density,  thus,  the  rate  of  rod  loss  is  reduced  and  the 
penetration  process  is  assisted.  In  the  second  zone,  the  apparent  density 
of  the  semi-infinite  material  is  greater  than  its  normal  density,  thus, 
increasing  the  rate  of  rod  loss  over  that  which  is  normally  expected  for 
a given  material.  One  reason  for  this  increase  in  apparent  density  is 
that  some  of  the  steel  plate  was  pushed  into  the  region  which  was 
initially  occupied  by  the  semi-infinite  material. 

The  thickness  of  material  in  the  second  zone,  necessary  for 
rod-target  interaction  to  be  completed,  depends  on  the  density  of  the 
semi-infinite  material  and  the  impact  velocity.  For  an  impact  velocity 
of  7 km/s,  rod-target  interaction  ceases  after  a penetration  of  about 
30. S mm  (1.2  in.)  into  a vacuum,  21.6  mm  (0.85  in.)  into  air,  and  17.5  mm 
(0.69  in.)  into  aluminum.  However,  for  an  impact  velocity  of  3 km/s, 
rod-target  interaction  ceases  after  about  5.1  mm  (0.20  in.)  of 
penetration  into  aluminum. 
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TABLE  V 


SHOCK  PROPERTIES  FOR  COPPER  ROD  - STEEL  PLATE  IMPACT 


Impact  Velocity 

Impact  Pressure 

Particle  Velocity 

Shock  Velocity 

km/s 

Mbar 

km/s 

km/s 

Copper  Steel 

3 

0.8 

1.57 

6.26 

6.34 

5 

1.7 

2.71 

8.02 

8.00 

7 

2.7 

3.67 

9.45 

9.30 

TABLE  VI 


SOU?4D  SPEED  AND  DURATION  OF  CONSTANT  PRESSURE 
Inpact  Pressure  Sound  Speed  t 


Figure  2 (U).  The  velocity  field,  at  7.345  usee  (computational 
cycle  #300)  after  the  Impact  of  a copper  rod 
on  a laminate  target.  Is  shown  with  an  outline 
of  the  rod-target  configuration. 
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particle:  velocity.  CM/mCROSEC 

Figurt  5 (U).  Shock  Velocity  - P«rt1cle  Velocity  Curves  for  Verlous  Meterlels 
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Interface  of  the  plate  and 
se«i-1nf1nUe  Material 


rm 


rod- target 


Figure  21  (0).  The  Motion,  along  the  axis  of  syametry  of  a rod-target 

configuration,  of  three  tracer  particles  is  shoMn  for  Case  T7 
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